Observing the Formation of Long-Range Order during Bose-Einstein Condensation 
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We have experimentally investigated the formation of off-diagonal long-range order in a gas of 
ultracold atoms. A magnetically trapped atomic cloud prepared in a highly nonequilibrium state 
thermalizes and thereby crosses the Bose-Einstein condensation phase transition. The evolution of 
phase coherence between different regions of the sample is constantly monitored and information 
on the spatial first-order correlation function is obtained. We observe the growth of the spatial 
coherence and the formation of long-range order in real time and compare it to the growth of 
the atomic density. Moreover, we study the evolution of the momentum distribution during the 
nonequilibrium formation of the condensate. 

PACS numbers: 03.75.Kk, 03.75.Pp, 05.70.Fh, 07.77.Gx 



When a gas of atoms is undergoing Bose-Einstein con- 
densation a macroscopic number of particles start to oc- 
cupy the same quantum mechanical state — it seems like 
the randomly colliding atoms are suddenly forced into a 
lock-step motion. The understanding of this process in 
which phase coherence spreads over the whole gaseous 
cloud has intrigued physicists long before Bose-Einstein 
condensation has been demonstrated [HQ. In particular, 
the question when the characteristic long-range phase co- 
herence is established is the key point for understanding 
the condensation process. The trajectory into the state 
of a Bose-Einstein condensate [H, 0, 01, 0] is much more 
intricate than its equilibrium properties. 

The transition to a superfiuid or a superconducting 
quantum phase is a remarkable process in which the prop- 
erties of the system undergo a fundamental change. The 
conceptual link between quantum phases in various sys- 
tems is the off-diagonal long-range order in the density 
matrix which describes phase correlations over macro- 
scopic distances [1, The first order spatial correla- 
tion function G^{r,r') = (*t( r )*( r ')) Q 

quantifies the 

characteristic length scale over which phase correlations 
exist. Here * (^ T ) denotes the annihilation (creation) 
operator of the atomic fields. To experimentally study 
how long-range order is established in space and time, 
real time access to this process is required. Yet, due to 
the short relaxation times and the strong coupling to the 
environment this seems troublesome in condensed matter 
samples, such as liquid helium or superconducting mate- 
rials. 

In a trapped atomic Bose gas the situation is distinctly 
different. It forms an almost closed system with neg- 
ligible coupling to the environment and the relaxation 
time scales are experimentally accessible. Previous ex- 
perimental studies of the Bose-Einstein condensate for- 
mation have focused on the increasing particle density 
[H B El and the influence of excitations on the momen- 
tum spread in elongated traps [13, E3 ■ The off-diagonal 
long-range order has been studied only at thermal equi- 
librium 12 1 . 



Here we present an experimental study of the evolu- 
tion of off-diagonal long-range order during the forma- 
tion of a Bose-Einstein condensate out of a nonequilib- 
rium situation. We start from a Bose gas above the phase 
transition temperature and suddenly quench the gas into 



a strongly nonequilibrium state [8|, |9|, llfj, [llj. Subse- 
quently, the gas can be regarded as a closed system which 
evolves into a Bose-Einstein condensed phase. The off- 
diagonal long-range order is measured by studying the 
interference pattern of atomic matter waves originating 
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FIG. 1: Formation of long-range order. Shown is the visi- 
bility of a matter wave interference pattern originating from 
two regions separated by Az inside a trapped cloud. This is 
a measure of the first-order spatial correlation function of the 
atomic cloud. By shock cooling the gas is prepared in a highly 
nonequilibrium state at time t = and then relaxes towards 
thermal equilibrium. Initially (t — 75 ms) the correlations are 
short ranged and thermal-like. The onset of Bose-Einstein 
condensation is marked by the appearance of long-range or- 
der extending over all measured distances (i = 275 ms). The 
largest measured distance corresponds to half of the Thomas- 
Fermi diameter of the final condensate. The data is taken 
for constant separations versus time. The substructure of the 
curves reflects the scatter of the initial conditions for each sep- 
aration. The resulting variance in the onset of condensation 
is about 30 ms. 
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from two different locations in the atom trap [l2j. The 
visibility of this interference pattern measures the phase 
coherence between the two regions and its temporal evo- 
lution uncovers the formation of long-range order. We 
continuously monitor the interference pattern during the 
formation using a single atom counter with high temporal 
resolution [l]| . Varying the vertical distance Az between 
the two locations allows us to experimentally map out the 
evolution of the phase coherence in the trapped Bose gas 
(see Fig.HJ. 

In the initial nonequilibrium state just after the quench 
we observe purely thermal-like short-range correlations 
which decay on a characteristic length scale given by 
the thermal de Broglie wavelength XdB ~ 0.4 pm. After 
the first 200 ms during which little changes are detected, 
the length scale over which phase correlations exist ex- 
pands rapidly and the long-range order of a Bose-Einstein 
condensate is established after a further 100 ms. Subse- 
quently, a stage of gradual condensate growth towards 
equilibrium is observed. In the Bose-Einstein condensed 
phase the system does — apart from its size — not exhibit 
a characteristic length scale 

In the theoretical investigation of the nonequilibrium 
formation of a Bose-Einstein condensate one distin- 
guishes qualitatively between a kinetic evolution and a 
coherent evolution 0, [l^, [lj|- During the kinetic 
stage of the evolution, the occupation numbers of the 
low-energetic states grow. It is governed by elastic colli- 
sions and the characteristic time scale is set by the col- 
lision time t co i = (navx)~ 1 with n being the peak den- 
sity of the gas, a the elastic collision cross section, and 
vt the average thermal velocity. The coherent stage — 
which plays a significant role only in the regime of very 
large density or scattering length — describes the merging 
of quasicondensates, patched regions of locally constant 
phase, into a full Bose-Einstein condensate. The presence 
of quasicondensates could be revealed by the momentum 
distribution of the gas below the critical temperature. 

Our experimental setup and the procedure for Bose- 
Einstein condensation have been described previously 
171 ]. For the experiments reported in this Letter, we 
start by preparing a thermal cloud of atoms in the 
\F = l,mp = —1) hypcrfme ground state of 87 Rb 
in a harmonic magnetic trapping potential. The trap- 
ping frequencies of the magnetic trap are (ui x , w y , uj z ) = 
2ir x (39, 7, 29) Hz, where z denotes the vertical axis. The 
temperature of the atom cloud T = 240 nK is slightly 
above the transition temperature T c — 215 nK for Bose- 
Einstein condensation for the given number of atoms 
1.3 x 10 7 . This results in a collision time of t co \ k 30 ms. 
Temperature and atom number are measured by absorp- 
tion imaging. The atom cloud is then brought into a 
strong nonequilibrium situation by rapidly lowering the 
trap depth to 620 nK and removing the high-energy tail of 
the Maxwell-Boltzmann distribution. Within the 100 ms 
of this "shock cooling" we remove 30% of the atoms. Sub- 



sequently, the cloud relaxes from its highly nonequilib- 
rium state and within a few hundred milliseconds 3% of 
the atoms form a Bose-Einstein condensate. During the 
relaxation process particle number and total energy are 
conserved with minimal disturbance due to the detection 
process. Moreover, we do not observe oscillations of the 
condensate with respect to the detector as a consequence 
of the shock cooling. 

We detect the evolution of both the density and the 
long-range order of the cloud simultaneously and in real 
time using radio frequency output coupling [l8| and sin- 
gle atom counting 1J, T3| • For output coupling we apply 
a weak monochromatic microwave frequency field which 
spin-flips atoms from the magnetically trapped state into 
the untrapped state \F = 2,rriF — 0). The untrapped 
atoms form a downwards propagating atomic beam. The 
output coupling region is defined by a surface of con- 
stant magnetic field and can be approximated by a hori- 
zontal plane within the atomic cloud [18]. Applying two 
microwave fields with different frequencies realizes out- 
put coupling from two vertically separated surfaces of 
constant magnetic field [l2| which are chosen symmet- 
ric about the center of the cloud. The two overlapping 
atomic beams interfere with each other. The mean flux 
reflects the local density of atoms at the position of out- 
put coupling within the transverse momentum interval 
measured by our detector. The visibility of the interfer- 
ence pattern reflects the phase coherence. The interfer- 
ence pattern is detected in time with single atom reso- 
lution using an ultrahigh finesse optical cavity, located 
36 mm below the magnetically trapped atoms [13]. The 
presence of an atom inside the cavity leads to a decrease 
in the transmission of a probe laser beam resonant with 
the empty cavity [20j | . Our system permits a precise mea- 
surement of the time-dependent atom flux. 

The mean atom flux and the visibility of the atomic 
interference pattern are determined in time bins of 50 ms 
length. In Fig. [2] three different situations of the con- 
densate formation are shown. For the cloud immedi- 
ately after shock cooling [Fig.^a)], the visibility is zero 
and the mean atom flux is low. As the condensation 
process develops [Fig.^b)], the atom flux increases and 
interference arises. Both continue to grow up to a fi- 
nal value [Fig.djc)] determined by the condensate frac- 
tion and detector function. The black lines are a fit of 
f(t) = A[l + V sin {2-wvt + <f>)] to the data. The fre- 
quency v is the difference of the two microwave frequen- 
cies used for output coupling. The phase 4> is determined 
by the relative phase of the two microwave fields which 
is locked to the experimental cycle. The output cou- 
pling of atoms with a rate « 5 x 10 4 atoms per second 
is essentially not perturbing the density dynamics of the 
condensate formation. In particular, it has no detectable 
influence on the final condensate fraction. 

The observation of long-range order even at a very 
small condensate fraction poses a severe technical chal- 
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FIG. 2: Buildup of density and off-diagonal long-range order 
during condensate formation. We plot the histogram of the 
atom arrival times modulo 1ms for v = 4 kHz. The data are 
summed over 18 repetitions of the experiment. The black 
lines are fits to the data to extract the mean atom flux A and 
the visibility V of the interference pattern (see text). 



lenge. Our detection scheme has specific properties which 
render the present measurement possible. The output 
coupling probability for a condensed atom is a factor of 4 
larger than for a thermal atom due to their different den- 
sity distributions for our equilibrium parameters. More- 
over, we observe unequal detection efficiencies for con- 
densed and noncondensed atoms which we estimate to be 
23% and 1%, respectively. They are mainly determined 
by the size of the atomic beam exceeding the active area 
of the detector. For thermal atoms the transverse veloc- 
ity spread is larger which reduces the probability for an 
atom to be detected in the cavity 



17j . Because of these 



effects the visibility of the interference pattern greatly ex- 
ceeds the condensate fraction and reaches a maximum of 
70% for Az > XdB- We observe a reduction of the visibil- 
ity towards large values of Az (Fig.[l|. This is influenced 
by the following effects: Output coupling symmetrically 
off center from the Thomas-Fermi profile leads to a re- 
duction of 7% at Az — 9 pm. Furthermore, the spatial 
modes of the two interfering atom laser beams depend 
on their divergence and thus on the position of output 
coupling [21| . We calculate this effect to be 10% for the 
largest Az. An additional contribution is due to the ge- 
ometry of our trapping potential resulting in a tilt of the 
long axis of the condensate by 4° with respect to the 
horizontal. Moreover, the residual uncertainty in deter- 
mining the arrival times of the atoms (4 us) limits the 
attainable contrast for the largest slit separations in our 
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FIG. 3: Growth curves for atom flux and visibility of the 
interference pattern during the formation of a Bose-Einstein 
condensate. The separation Az between the two output cou- 
pling regions is 1.9 um. Before the shock cooling stage there is 
a low flux of thermal atoms and no interference. The growth 
in atom flux is significantly slower than the growth of the in- 
terference pattern. The data are averaged over 20 repetitions 
and the error bars indicate the statistical error of these mea- 
surements which are within a factor of 2 of the theoretical 
shot-noise limit. 



experiment to 94%. 

In Fig.[3]we show the growth of the density reflected in 
the mean atom flux and of the visibility for a given sep- 
aration of the output coupling regions. This corresponds 
to a section of constant Az = 1.9 pm in Fig.Q] We have 
analyzed both the delay of the formation with respect to 
the shock cooling stage and the speed of the formation 
for density and off-diagonal long-range order. We fit a 
function g(t) = A 2 + (Ai - A 2 )/[l + (t/r) p ] to the data 
to quantify the growth, where r denotes the time after 
which 50% of the increase in flux or visibility are reached. 

The duration of the condensate formation can be quan- 
tified by the time needed for an increase of both the flux 
and the visibility from 10% to 90% of the total increase. 
We find the duration to be approximately independent 
of the separation Az of the two output coupling regions. 
For the flux this time is 421 ± 47 ms whereas the visibil- 
ity of the interference pattern grows faster in a time of 
267 ± 48 ms. From the data we find that the time r after 
which 50% of the total atom flux is reached is almost con- 
stant across the investigated inner region of the trapped 
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FIG. 4: Temporal evolution of the momentum distribution of 
the detected atoms. The distribution is obtained by Fourier 
transformation of the visibility data. The thermal distribu- 
tion is calculated for the parameters before shock cooling in- 
cluding a normalization factor. 



cloud, whereas for the visibility it increases with larger 
slit separation. The coherent region grows with a veloc- 
ity of « O.lmm/s, approximately a factor of 5 slower 
than the speed of sound at the peak density of the ther- 
mal cloud. The speed of sound imposes a natural speed 
limit for the expansion of the coherent region. While the 
onset of the growth (defined by the 10 % points) starts 
39 ± 54 ms earlier for the flux than for the coherence, the 
visibility saturates (90% points) about 116 ± 34 ms ear- 
lier than the density. This applies to all separations Az 
larger than the de Broglie wavelength. 

Moreover we find that r decreases with increasing size 
of the final condensate. To obtain larger condensates in 
equilibrium, the initial particle number or temperature 
or the fraction of atoms removed during shock cooling 
was changed. As the limiting cases we find for a final 
condensate number of 2.5 x 10 5 atoms r = (282, 320) ms 
and for 8.9 x 10 5 atoms r = (57, 114) ms, where the first 
number in the bracket refers to the visibility and the 
latter to the flux. The numbers are averaged over all 
slit separations. This decrease of the formation time for 
larger condensates is proportional to the increase of the 
elastic collision rate due to the higher final density. 

It is instructive to study the evolution of the vertical 
momentum distribution n(k) during the thermalization 
process. We infer the distribution n(k) of the detected 
atoms from the Fourier transform of the visibility data 
of Fig. [T] multiplied by the average flux (2]|. The result 
is shown in Fig.[4j For the initial state directly after the 
quench, the momentum distribution is well approximated 
by a thermal distribution (dashed line) . The thermaliza- 
tion process after the shock cooling leads to the formation 
and the growth of a low-momentum peak which in equi- 
librium corresponds to the condensate. For our parame- 
ters in equilibrium we expect a true condensate without 



strong residual phase fluctuations [23|]. In Fig.|4j we fo- 
cus on the beginning of the condensate formation. After 
225 ms of evolution, we observe a broader low-momentum 
peak (width in fc ~ 0.8 urn -1 ) than for the final distribu- 
tion. For evolution times of 325 ms and larger, the width 
of the low-momentum peak is limited by our measure- 
ment resolution of 27r/(16um), which takes into account 
that the correlation function is symmetric with respect 
to Az = 0. 

In summary, we have experimentally studied the dy- 
namics of Bose-Einstein condensation and more specif- 
ically the formation of long-range order using an inter- 
ference measurement. We have quantitatively compared 
the evolution of the flux and of the visibility of the inter- 
ference pattern. Finally, we have presented our data in 
terms of the atom momentum distribution. 
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